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From the theory of quantum LC circuits with discrete charge, and semiclassical considerations, 

we obtain approximate energy eigenvalues, depending on the parameter q^/h. Next, we include 

qq ' electrical resistance for the quantum RLC circuit, obtaining a relation that strongly reminds us of 

f*^ ' the Landauer formula. 
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PACS numbers: 73.21.-b, 73.23.-b, 73.63.-b 



I. INTRODUCTION 



In a series of articles Li and Chen 0, J2| and us 0, |j, [3 m 0] , have developed a theory of quantum electrical systems, 

based on a treating such systems as quantum LC circuits; that is, electrical systems described by two fenomenological 

parameters: an inductance L, and a capacitance C '. As it is known, when the transport dimension becomes comparable 

with the charge carrier coherence length, one must take into account not only the quantum mechanical properties of 

' ■ the electron system, but also the discrete nature of electric charge, which leads to the concept of quantum LC circuit 

qj \ with discrete charge. 

In this work, we devote section [TT] to a brief discussion on some basic facts about quantum circuits with discrete 
charge, to set up the background for this work. In section IIIII we discuss the semiclassical approach to the study of 
quantum circuits, computing the approximate energy eigenvalues using the Bohr-Sommcrfeld quantization rules, for 
two different energy regimes. The approximate energies are expanded in a power series on the dimensionless parameter 
y/L/C q^/h. In section IIVI we generalize the equations of motion for the RLC circuit, the resulting equations are 
'"O ' equivalent to the damped simnple pendulum. Finally, we find a set of stable solutions, characterized by constant 
charge and flux, finding a simple relation between the RLC circuit parameters that bears a striking resemblance to 
O ■ the famous Landauer formula [3 Q • 



II. SUMMARY OF QUANTUM CIRCUITS 



> 

qq ' In their pioneering article, Li and Chen [l|, [2| consider a quantum LC circuit, described by a wavefunction ^f(q) in 

^f , the (continuum) charge (q) representation, 

g; H*(q)=- — — *(q) + V(q)*(q). (1) 

00 ' 

In this description, the charge (q ) and flux ($) operators, satisfy the usual commutation rules, [q, $] = ih, and may 
;> be represented by the operators q = q, and $ = —ihd/dq. The charge operator posseses continuum eigenvalues (go), 

with delta function eigenstates, just as the single particle quantum mechanical operators q and p. 



qy qo (q) = q y qo {q)=q Q 6(q-q a ). (2) 

Li and Chen introduce discrete charge eigenstates, whith eigenvalues equal to integer multiples of the quantum of 
charge (the electron charge, g e ), g„ = nq ei 

g*„(g) =ng e *„(g). (3) 
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In the charge representation, the Schrodinger equation, H'fy(q) — E^(q), becomes a difference equation, 

2 [*(g + q e ) + tf (« - q e ) - 2tt(g)] + V(q)*(q) = EM(q), (4) 



h 2 



2Lql 

corresponding to a discrete hamiltonian H = Hq + V(q), defined below. It is convenient to define [l| the discrete 
charge-shift operators Q and Q', and a discrete flux operator <t, 



0*(<z) = 


= *(<7 + <7e) 


Q+*(g) = 


= *(?-g e ) 


# o = 


-2^ (Q + Q '- 


$ = 


= (B/2w.)(<5-<5') 
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It is found that the discrete flux operator $ satisfies modified commutation relationships with q and H below, due 
to the discrete nature of electric charge, 



q, $ 


= ih(l + ^H 


H ,$ 


= 


H ,q 


= ih<&. 



In the charge representation, the operator <t posesses "plane waves" eigenstates, that is solutions of the equation 
<S>iS>(q) — X^/(q), given by ^(q) — exp(i<fiq/h), with eigenvalue A = sm(q e <p/h)/(q e /h), in which is a continuum 
quantum number, which we call the pseudo flux |5|. With this one may define the " pseudo-flux" representation, in 
which one deals with " wavefunctions" 4" (</>), the charge operator becomes a differential operator with respect to <j>, 
q = ihd/d(f>, and hence the hamiltonian operator becomes a differential operator also, just as in the continuum charge 
case, 

h 2 d 2 2h 2 9/ , s 
H = -2CW + W^ {qAm - (6) 

Observe that the physical flux is the operator $ = (h/q e ) sm(q e (j)/h) ^ <f>, then, the operator 4> will be called pseudo- 
flux from now on. Note also that the physical current operator is obtained from the canonically invariant definition 
I = — i[H, q]/h, which gives I = <t/L, which is why we consider $ to be the physical flux. 

III. SEMICLASSICAL STUDY OF LC CIRCUIT 

The electrical engineer has learned to love the simplified description of a system provided by the circuit approxi- 
mation, when compared with the more complete, but also more complex field description. On the other hand, nature 
is quantum, we say; however, we describe the behaviour of electrons in modern circuits using the same basic laws 
(Kirchhoff laws) as in a classical circuit. There are difficulties on sight for this state of things, since we are now prob- 
ing nature at very low temperatures, with very pure materials, very tiny currents and strong magnetic fields. Many 
examples show that things are about to change for the engineer, we have seen flux quantization on superconductors, 
conductance oscillations, quantum hall effects (integer and fractional), persistent currents and so on [alSl- 

It would be very useful to find out to what extent a circuit-like description could be of use for the very small 
electronic circuits of tomorrow, and what may be retained from such description; a description that has the value of 
simplicity and familiaruty, which we [3|, |4j, [5|, |g, [7[ have been pursuing for some time. Undoubtely, such a descritpion 
will be able to answer some questionsonly, leaving others unanswered. For example, one area in which the "quantum 
LC circuit" may give valid results is in the calculation of energy spectra. We have done this for the open electron 
resonator 0, |T(| ; but one may question that the calculation was not " simple enough" , since it still required to solve 
the Schrodinger equation. 



Now, we propose to go one step further in our simplification, by proposing to use a " semiclassical" approach to 
the study of quantum circuits. Note that the word semiclassical will be used in a loose sense, as we shall see. We 
start from our model hamiltonian in the pseudo-flux representation (equation [6] above), and treat all the operators as 
classical variables; therefore we shall drop the 'hat' from the hamiltonian (H — > H) and the variables (q — ► q, (j> — > 4>), 

a 2 2h 2 
H = ^ + — S m\q e <f>/2h). (7) 

The hamiltonian equations corresponding to eq. ([7]) coincide with the Heisenberg equations for the quantum hamil- 
tonian above, but with all the 'hats' suppressed, i.e. 

dH q , 

"da- = C = -+ 

_ = _ an^/h) = q, 

Observe that H is not a true classical hamiltonian; for example, formally letting h — > yields nonsense, unless one 
takes first the limit q e — > 0, which would be nonsense for us too. Let us define 0o = h/q e , noting that the flux quantum 
tpo = h/q e = 27T(po, so that the hamiltonian and the corresponding equations look classical, since we eliminated the 
explicit reference to Plank's constant, 

H = lU^sin^^o). 

V C 

q = —sm(4>/(/)o). 

The classical equations above (q and <j) are c- numbers), may be written as a single second order equation, for 
the dimensionless variable 9 = <j)/4>o, which coincides with the equation for the simple pendulum, with frequency 
f2o = 1/y/LC, = — £1q sin(^) , which has been extensively studied form many years, becoming a standard textbook 
exmpla in classical mechanics |ll| . We are interested in this classical equation, as a means to obtain the eigenstates 
of the quantum hamiltonian eq. [BJ the method is the Bohr-Sommerfeld quantization rules (WKB approximation) of 
the old quantum theory, i.e., § qdtfi = (n + j)h (7 is a constant, we take 7 = 1/2). 

Observe that the constant (J)q defines a scale of energy for our quantum LC system, namely <pf)/ L, therefore,, there 
exist two energy regimes, the low energy regime, defined by E < 2</>§L, and the high energy regime, defined by 
E>2(j) 2 /L. 

A. Low energy states: E < 20QL 
We consider here the calculation of the energy spectrum of the quantum LC system, in the 

q = CJ> = Ccj> e, 
E = Ccf> 2 I y + 2Q;j sin 2 (0/2) J . 

Let 9 = 9o be the maximum amplitude, defined by sin 2 (# /2) = E / (2C (J) 2 Qq) = E/(2cf%/L), then 



(9 = ±2ft y / sin 2 (V2)-sin 2 (6y2). (8) 

To find the (quantum) energy eigenstates we need the phase-space integral /, 

I = j> qd4> = C<j>lj)9d9 = 2 J j(j)l I ^/sin 2 (0 o /2) - sin 2 (6»/2) dS (9) 



which may be expressed as an elliptic integral of the second kind, expressed in the standard form by the transformation 
sin(£) = sm(#/2)/sin(#o/2), defining further k — sin(6*o/2), we get 



I = mJ^lk* f ' -1 
v l Jo u 



y/l-k 2 sm 2 (£) 



This integral is given in the tables [T4], page 162 (2.584, #6), in terms of the complete elliptic integrals, K(fc) 
F(tt/2, k) and E(fc) = E(tt/2, k), 



I = I6\l j^ (E(k) - (1 - k 2 )K(k)) 



(10) 



To obtain approximate expressions for the energy eigenvalues, we use the series expansion for K(fc) and E(fc) (Grad- 
shteyn [12j |. pp. 905), valid for k « 1, we get 

/ « A^cgk 2 (l + £ 
Impose now the quantization rules of the old quantum theory, I = (n + ~f)h, to obtain the energy from 



iv^M^U+y) = („ + ~.i/,. 



Assuming k « 1 we have, as a first approximation, 



(11) 
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which gives the usual result for the energy, E® = (n + j)hQ,Q. Now compute the corrected k n , 



(n + j)h [L_(_l{n + i)h jL 



kn 4%$ V C y 8 4tt$ V C 
The energy levels become (using 7 = 1/2, as indicated previously), up to second order 



E„ 



hfl Q (n + j) 1 



hQ (n + 7) 



1 (n + j)h [L 



8 4tt02 V C 
n(n + "f) q 2 L 



h V C 



(12) 



This result becomes E n = (n + l/2)M2o) m the continuum charge approximation q e — > 0. Notice also that, on the 
next approcimation, the result depends on the ratio of the Landauer conductance q 2 /h, and LC conductance y/C/L, 
an interesting result, as we shall see later on. 

B. High energy states: E > 2<f>o/L, or C-design case 

The so-called C-design case corresponds to the 'large capacity' case, in which the electrostatic part dominates. We 
have seen that the 'velocity' <fi may be written as 



V c 
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Consider the case in which the energy E » (f)^/2L, i.e., in which the electrostatic energy dominates, the so-called 
C-design. It is convenient to define the parameter A = (2</)q/L)/E, to write the expression for the electric charge q, 



q = ±V2CE 



1-Asin 2 (-^) 



1/2 



(13) 



Observe that in the present case, the pseudo flux variable cj> grows without bounds, but the system is periodic with 
period A0 = 2tt<J)q. To compute the action integral /, it is convenient to change variables to a — (j)/2<fio, 



I = Iqd(f> = V2CE2<f>o J [1- Asin 2 (a)] 1/2 da = 2(J) V2CEF,(tt, VA), 



(14) 



in which we notice that and that E(tt, \/A) is an elliptic integral ([12| equation 8.121.4, pag. 906), that E(7r, \/A) 
2E(tt/2, VX) = 2E(v / A), and also 2(/) Q V2~CE = (i4^/VX)y/U/L, hence, we have 



Now, use the series expansion of E(v / A) (uMi equation 8.113, pag. 905), and replace (f> — h/2nq e 



(15) 



I = 
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Now, impose / = nh (use 7 = for the Maslov index), and compute the first approximation Ao, 



(16) 
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(17) 



The condition A < 1 becomes n > y/ 'C '/ 'L / '(irq^ / 'h) , for the quantum number n. This is interesting, since y/C/L 
is the conductance associated to the LC circuit, and q 2 /h is the conductance quantum, appearing in the Landauer 
formula, this ratio appears as an expansion parameter. Define the dimensionless ratio fj, = A/Ao, then the following 
relation is equivalent to the quantization condition / = nh, 



VJi- 



1 22 6 8 2 16 2 



(18) 



Now, assume the expansion /_* = 1 + aAo + b\$ + cAg + • • •, and insert into the equation above (after squaring both 
sides); the result for A becomes 

A = Ao (1 - 2A + 29A 2 ,/32 + • • •) , 

A" 1 = A^ 1 (1 + 2A + 99A^/32 + • • •) 

then the approximate energy eigenstates are 



= 24>l _ 24% 40g 990 2 Aq 
LX LX L 16L 

= (ng e ) 2 m 99 f^~ 

2C L 32 V L J (nq e 



C 



(19) 
(20) 



IV. ELECTRICAL RESISTANCE AND CIRCUIT LAWS 



The problem of electrical resistance at the mesoscopic level is a very subtle one. The problem has been studied 
by Landauer and others [8j, who obtained a relationship between conductance and the transmission coefficient in a 
one-dimensional scattering experiment. The role of the contact resistance, as well as that of the different physical 
regimes in which one may observe a system, complicate the description of a system. 



A. The generalization 

Here we propose to study a generalization of our equations to include electrical resistance, in a manner similar to 
what happens in a classical circuit. To carry out this generalization, we may proceed by replacing the physical current 
/ occurring in a mesoscopic circuit by the value X, and adding a 'resistance' term, in the same way as it is done in a 
classical circuit, i.e. by adding a potential drop AV = RX = R(j>/L; then the equations become 

-1-BX = 1 (21) 

q = —sm(c/)/(j> ). 

In this way, we have a generalization that preserves the form of Faraday's law of induction, e = —d(j)/dt, and Ohm's 
law, AV = EX; which is rather nice for engineering applications. Now, we must decide which value should we assign 
to the X appearing on equation 1211 We note that there are two possible choices for it, namely, we may identify 

• X = q, or 

• X = 4>/L. 

Let us consider first the (possible) identification X = q. In this case, the system is described by 

-4>-R^M4>/M = § (22) 

q = —sin((j)/(j) ). 

If we study the equilibrium solutions of eq. 1221 we see that q = implies cf> = rnrtpo/^, and that = 0, implies that 
the charge q = 0. 

Consider now the second choice, our proposal, in which X = <f>/L. The quantity X may be called a pseudo-current, 
since <f> is already called pseudo-flux by us; notice then that X =/= q, which is ok, since the true current is I = q. Our 
new circuit equations become 

-t-4 - j m 

4 = -j^sm(<j)/(l) Q ). 

To make our generalization more compelling, we prove that it may be obtained from a lagrangian formulation. To 
see this, we write the equation of motion as a single, second order equation, 

L4> + R<j)+^ sin((/)/0o) = 0. (24) 

This equation may be obtained from the Lagrangian description, with the use of a " dissipation function" , as 
described by Landau [llj. Let C be the lagrangian, and T be the dissipation function, then 

£ = ^-^™ 2 »/2« (25) 

r - *£ (2 6 ) 

— ( — )- — -- -— (27) 

dt K d<j> d<j> dj> 

The last equation then is exactly our proposed equation. This is a good generalization, in the sense that the energy 
E, defines in the usual way, 



E=^f- + ^sm 2 (0/20 o ), (28) 



decreases, or stays constant (if = 0), since 



L 



- = -W. (29) 



The new equations have also the correct classical limits, obtained for q e — > and small pseudo-flux ( <fi << 0o), 
they also have more interesting structure; for example, the equations posess nontrivial solutions with constant pseudo 
flux and constant charge, to be discussed below. 

V. RELATION TO LANDAUER FORMULA 

As it was indicated above, the system of equations [23] posseses constant flux and charge solutions, which we find 
by assuming q — and <j) — 0, we obtain (j) n /4>o = nw, for integer n, and we get the condition —R<p n /L — q n /C, 
i.e., q n — —RC(j) n /L — —nnRCfio/L. The solutions having even n are stable, and those with odd n are unstable, as 
it is shown by numerical calculations: the solutions for odd n "decay" onto the n — 1, i.e., dissipating energy in the 
process, while 0^0. Hence, from now on we keep only the stable, even n solutions (replacing n by 2n). We have 
shown then, that there exist stable equilibrium solutions of the equations of motion [2"31 which may be expressed as 

<f) = 2mr(f>Q 
q = — 2nq* 

TT(j>oRC 

q = — — 

These solutions show that the pseudo flux <j) is an integer multiple 27r0o, and that the charge is an even integer 
multiple of the charge q* — ^qRCJL. The physical, equilibrium charge in the system is q; therefore, this charge q 
must be an integer multiple of the unit of electronic charge (q e ), i-e., it should be written as q = m(— q e ), in which m 
is another integer number, not necessarily equal to n, as we shall see below. We have 



2nRC-K(, 



-mq e = , 

then we obtain 

RC m q e m 2q\ m 

7 — 7\ ] rt t Landauer- 

L In 7r0o In h n 

The quantity Ghandauer — Qe/h — 3.875 x f0~ 5 Q,^ 1 is the Landauer conductance. This is a remarkable result, 
obtained by a simple generalization of the equations of motion to account for electrical resistante. Notice also that the 
quantity C/L is actually the product of the C— conductance and an L— conductance, in other words, recalling that 
the capacitive and inductive conductances at frequency flo, Gc — QqC, and Gl = 1/QqL (at frequency ui) satisfy 
GcGl — C/L, we may write, or that Glc = \JGJL may be viewed as the conductance of thew LC circuit. 

G e ff = RG LC — — = — G Landauer- (30) 

n h n 

We may say that the circuit should show an effective conductance G e ff, which is a rational multiple of the Landauer 
conductance Gl- We have established a relationship between the parameters of a mesoscopic system, that bears a 
strong resemblance to the Landauer formula, this result comes about because we have insisted on charge quantization. 
This is analogous to what happens in the old quantum theory, in which the angular momentum quantization, applied 
to the Bohr atom actually selects the energies. In our case, our conditions selects the cases in which we have both 
flux and charge quantization, resulting in equation 1301 

The previous relation may be generalized slightly if one considers the influence of an external, constant electromotive 
source eq, in this case, the circuit equations become 



^ + ^ + R l ( 31 ) 

— sin(0/0 o )- 



Now, let us impose the condition q = 0, and <f> — 0, then = rvK(j)Q, as before, and the charge becomes 

R<f>o~ 



C 



£o — nit- 



L 



(32) 



Now, as in the previous paragraphs, we consider only the stable (even n) solutions, i.e., replace n — > 2n, and impose 
the condition q = —mq e , for integer m, then we obtain the relation 

mGzandauer = — CVo^O + n G e ff (33) 

Finally, if one considers the system as driven by an external, time-dependent perturbation, v(t), the equations of 
motion may now be written by simply adding the external source, as below 

* = -f-%+v(t) (34) 

q = —sm(c/)/(f> ). 

It it clearly seen it is not possible to obtain an equation fot the charge variable q alone, but that it may easily 
done for the pseudo flux variable <f>, the result has a simple and well known form, namely, that of the forced simple 
pendulum with dissipation. These equations describe may be used to describe many different physical systems, 
including mechanical, electrical, superconductors, and so on, and it will be considered elsewhere;for the variables of 
our semiclassical RLC circuit the equations are, which we write for the sake of completeness 

0+^ + ^sin(0/0 o ) = ^). (35) 

VI. FINAL REMARKS 

We started from a description of the quantum LC circuit, and introduced the semi classical approximation to 
compute the energy eigenvalues of the system. We computed the energy eigenvalues for two different energy regimes. 
First, we considered the low energy regime E < <Pq/L, and next, the high energy regime E > <\>\/L. In both cases we 
computed the energy as a power series in a small dimensionless parameter, which may be defined as y/ L/C (q^/h), 
i.e., the Landauer conductance times the impedance y/ 'L/C of the LC circuit. In this way, we have shown in a very 
simple way the important role of q^/h, and its relation to charge quantization on the LC circuit. 

Next, we generalized the equations of motion of the LC circuit to account for electrical resistance, which is equivalent 
to generalizing the usual Kirchoff circuit laws, to account for charge discretization and electrical resistance; the 
equations obtained arc equivalent to the simple pendulum with forcing. The sin(0/0o) term in the equations of 
motion generates a set of stable constant flux and charge solutions for the unforced mesoscopic RLC circuit, which, 
when charge quantization is imposed on the solutions, generates a condition to be satisfied by the parameters of 
the system 1301 This equation bears a striking resemblance to the Landauer formula and also to the quantization of 
conductance of the quantum Hall effect. This is a most amazing result, since it has been obtained from a generalization 
of the circuit equations, and simple considerations of charge quantization. 
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